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ON THE ANTIMATTER SIGNATURES OF THE 
COSMOLOGICAL DARK MATTER SUBHALOS 



Julien Lavalle 



Abstract. While the PAMELA collaboration has recently confirmed 
the cosmic ray positron excess, it is interesting to review the effects of 
dark matter (DM) subhalos on the predicted antimatter signals. We 
recall that, according to general subhalo properties as inferred from 
theoretical cosmology, and for DM with constant annihilation cross 
section, the enhancement cannot be > 20 for the antimatter yield. This 
bound is obviously different from that found for 7-rays. We also recall 
some predictions for supersymmetric benchmark models observable at 
the LHC and derived in the cosmological N-body framework, showing 
in the meantime the existing discrepancy between profiles derived from 
N-body experiments and the current observations of the Milky Way. 

The positron (e + ) cosmic ray (CR) fraction has recently drawn considerable 
attention due to the confir mation of a e + excess reported by the PAMELA col- 
laboration up to 270 GeV ( Adriani et al. . 2008). While the amplitude of this 
excess is unclear a ccording to the most recent predictions of the secondary e + 



flux at the Earth (jDelahave et all I2008L where we have emphasized the role of 
the e~'s), dark matter (DM) annihilation has been again proposed as a source of 
such e + 's, as it was already the case with the less clear HEAT excess ten years 
ago ( Barwick et al. . Il997t iBaltz fc Edsiol 19981 ). However, this primary contri- 
bution usually needs to be amplified to fit the data, and most of the authors 
have invoked some ad hoc boost factor coming from the presence of DM substruc- 
tures in the Milky Way (MW), which should indeed increase the average anni- 
hilation rate. Such subhalos are expected in the frame of hierarchical structure 
formation, and well resolved down to ~ 10 5 M Q in cosmological N-body simula- 
tions as performed in the A-CDM scheme. Here we shortly review the effects 
of DM clumps on the antimatter yield, focusing on e+'s and antiprotons (p's). 
We will not discuss the PAMELA e + fraction, since it is worth waiting for the 
release of the e^ spectra t hemselves. The results an d figures that we present 
here have been derived in ( Lavalle et al. . 2007 . 2008bl la). to which we refer the 
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reader for more details and references. Relevant scales for e + 's and p's. As 
the DM distribution is usually found to exhibit a steep spatial dependence in N- 
body simulations, the knowledge of the typical CR propagation scale is useful: it 
characterizes the spatial extension over which the observed signal at the Earth 
is integrated, and depends on both the propagation model and the CR species. 
In the GeV-TeV energy range, CRs diffuse on the Galactic magnetic turbulences, 
which ensure their confinement in the MW in a region which can be extended by 
a few kpcs above and below the Galactic disk. They interact through various pro- 
cesses with the interstellar medium (ISM) and/or the interstellar radiation field 
(ISRF). When propagating, p's mostly experience spallation in the disk, where 
the ISM is located, and convection outward the disk. These processes are efficient 
at low energies < 5 GeV, while energy losses are almost irrelevant; reacceleration 
can also be neglected. e + 's obey to the same spatial diffusion as p's, but lose 
very quickly their energy, mostly through scattering on the ISRF. Their propaga- 
tion is thus mainly set by the synchrotron and inverse Compton energy losses, of 
typical timescale r ~ 300 Myr, which occur in the whole diffusion zone. Other 
processes can be neglected without loss of accuracy above a few GeV. This al- 
lows to infer the characteristic propagation lengths for e + 's and p's: Xp — K \f^ 

and X 2 e+ = {(H)" 5 1 - (if-) 5 I' where y c ~ 12 km/s is the convection 

velocity, and E$ is the e + injected energy at the source, K(E) = Kq(E / Eq) s is 
the diffusion coefficient, with S ~ 0.7, Kq ~ 1.12 x 10 _2 kpc 2 /Myr being the nor- 
malization at Eq = 1 GeV. The propagation length for a p quickly increases with 
energy scaling like the diffusion coefficient oc E s . This is the opposite for e + 's, 
since the propagation length increases only as a e + loses its energy. This is illus- 
trated in the left panel of Fig. [TJ For e + 's, the characteristic propagation length is 
of order of a few kpcs, much lower than for p's, so the signal at the Earth comes 
from limited regions around, slightly enlarging as the detected energy decreases. 
Large integration volumes for e + 's (p's) will be found at low (high) energy only, 
even reaching regions close to the Galactic center (GC). Therefore, local fluctu- 
ations of the injection rate mostly affect the high (low) energies for cosmic e + 's 
(p's). DM distribution in the MW and subhalos. The DM annihilation rate 
is set by the squared DM density. In the frame of A-CDM cosmology, the DM 
distribution is found to be rather scale invariant and close to universal in N-body 
si mulations, when baryo ns are not considered. This is the seminal result obtained 
bv lNavarro et al. ( 19971 ) (NFW), and confirmed, though with some scatter in the 
parameters, by many other groups since then. These scale invariant profiles are 
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often given in the following spherical form: p(r) = po 

where r s is a scale radius, the index will refer here to the Sun location, such 
that ro = Rq — 8 kpc and po = P© = 0.3 GeV. cm -3 . A NFW profile corresponds 
to (a,/3, 7) = (1,3,1), and r s is usually found ~ 20 kpc for MW-like objects. 
Although it makes sense to constrain the DM distribution from a top-down ap- 
proach, it might sound risky to use such profiles, obtained without baryons, to 



On the antimatter signature of the cosmological dark matter subhalos 



3 



make predictions at the sub-galactic scale, where baryons are expected to play a 
major role, but this has been done by many authors (including ourselves). Indeed, 
some observational constraints exist, so do many co nsistent dynamical modelings 
of the baryon component of the MW. In particular, lEnglmaier fc Gerhard! (|2006h 
subtracted their baryon mass model to the velocity curves derived from CO and HI 
observations, and provided a set of data points on which one can directly constrain 
the DM contribution. In the right panel of Fig. [2], we show these data, on which we 
draw several density profiles. One is the NFW profile as employed for the MW, the 
others h ave been derive d from fits on a N-body MW-like galaxy of the HORIZON 
project ( Tevssierl l2002t) . These profiles clearly overestimate the central DM mass, 
so the DM annihilation rate close to the GC as well. While this will not affect 
that much the predictions for antimatter CRs, which are mostly set by the local 
density, this is instead expected to lower the 7-ray predictions. Now, we shortly 
turn to DM subhalos. For what concerns indirect detection of DM, the important 
features are their predicted spatial distribution, mass distribution and inner indi- 
vidual profile. There is no unique values predicted for these ingredients, so a way 
to estimate the theoretical uncertainties is to bracket each of them with the most 
extreme possibilities. As regards the spatial distribution, it is found antibiased 



with respect to the smooth profile in some N-body experiments (jDiemand et al 



20071 ). which means P c \{r)/p s (r) cx r, but we have also studied the case where 



subhalos track the smooth profile. The mass distribution is usually found close 
to M~{ aM ~ 2 \ as expected in the linear theory, and olm S [1-8,2] seems reason- 
able. The minimal mass is set by the free streaming scale of DM particles at 
the collapsing time. For weakly inter acting particles with m asses ~ 100 GeV, the 
minimal clump mass is ~ 10" 6 M Q (|Profumo et all l2006h . which will be taken 
as a minimum here. The mass function can be assumed to be independent of 
the spatial distribution, so that the subhalo number distribution will be given by 
dn/dM = N tot ■ dV c \{M)/dM ■ dV{r)/dV, where 2V tot = (M d )/M c i tot is the to- 
tal number of clumps in the MW. The inner subhalo profile is mainly featured 
by its inner logarithmic slope 7 and the so-called concentration parameter. The 
latter is defined by c v j r = r V i r /r_2, i-e. the ratio of the subhalo virial radius to 
the radius at which the logarithmic slope of the profile is -2, and is found to de- 
crease with the subhalo mass. Different models exist, usually fitted on N-body 
data, and we use two extreme configurations w hich encompas s the w ide range of 
possibilities: one optimistic e xtrap olated from iBullock et aD ( 2001 ). and a sec- 
ond minimal from lEke et al.1 ( 2001 ). referred to as B01 and ENS01. Finally, the 
main quantity at stake when computing the antimatter flux is the intrinsic sub- 
halo annihilation rate, which depends on its mass, concentration, and inner profile: 
£(M c i,Cvi r ) = J j d 3 x(p c \/ po) 2 , which has the dimension of a volume, and is nor- 
malized to the local rate. Note finally that clumps are not expected to be numerous 
in the central regions of the MW because of tidal disruption, and in any case not 
concentrated enough to dominate over the smooth annihilation rate there. Fluxes 
and boost factors. The antimatter CR flux originating from the smooth DM 
halo is given by <f> s {E) =Sj d 3 x {p s (x )/p a ) 2 - J E dE s G(x Ql E «- x, E s ) ■ dN/dE s , 
where S cx (av) (p / m x ) 2 depends on the DM particle properties (whose annihila- 
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tion cross section is assumed constant), Q(f <— i) is the Green function associated 
with the CR propagation, and dN/dEs is the injected CR spectrum at the source. 

The flux due to subhalos is given by #f (E) = Ci(M c i, c vil )g(x Q <- x t ) 
S (N^) (£) (Q (x©)) , where Q(f <— /) is the convolution of the Green function with 
the injected spectrum, and () means the average over the relevant subhalo distribu- 
tion (space or mass). The boost factor is given by B(E) = (l-/) 2 + 0*° t (i?)/^ s (S), 
where / is the fraction of the local density in forms of clumps. It does depend 
on the energy, because of the energy dependence of the propagation scale, which 
sets the integration volume. It is independent from the injected spectrum only 
for p's. It is useful to derive the local asymptotic limit of the boost factor, since 
this allows the computation of its maximal value. It is valid at short propagation 
scales, i.e. at high (low) energy for e+'s (p's): B Q ~ 1 + (iV^ *) (^)dV{x Q )/dV , 
where / is neglected. This expression does not depend on the CR species, and 
provides a good way to check numerical computations. We have calculated the 
boost factors and associated statistical variances for p's and e +/ s, by considering 
a flat spectrum for the former, and a 200 GeV monochromatic line for the latter. 
We have studied different subhalo configurations, and different propagation mod- 
els. The latter test mainly influences the variance, slightly the energy dependence 
of the central values, since it only affects the propagation scales. The former is 
instead responsible for a huge scatter in the single clump luminosity, but, inter- 
estingly, the central values still lie in a small range, since the smooth component 
is found to dominate the overall signal in most of cases. As shown in the middle 
and right panels of Fig. [TJ the maximal value obtained for both e + 's and p's is 
~ 20, which corresponds to the maximal clump configuration: ~ 10 16 clumps in 
the MW, of minimal mass of 10 _6 M©, steep mass distribution cx M c \~ 2 , spatially 
tracking the smooth NFW halo, with r~ 3 / 2 inner profiles and B01 concentrations. 
Such a configuration is theoretically very optimistic, so rather unlikely. Taking 
clumps of 10~ 9 M© would give an asymptotic value of 23© ~ 30, which is still very 
weak, and considering a nearby very massive clump is unlikely and observationally 

3/2 

constrained. The relative variance scales like ~ A d , as expected: the number 
of sources is greater in larger volumes. Moreover, at large Ad, the boost is close to 
1 and the variance is small in any case, since the smooth contribution of the GC 
comes into play and dominates the signal, further shrinking the variance accord- 
ingly. We finally went beyond this analytical study by directly using a N-body 
galaxy from the HORIZON project, in which we could study more complex effects, 
e.g. density fluctuations (beside clumps themselves) or departure from spherical 
symmetry. We have quantified these effects to be small, and shown that a mere 
spherical smooth halo provides a good estimate of the antimatter fluxes. We show 
the predictions that we have derived for supersymmetric benchmark models ob- 
servable at the LHC in Fig. [2] Why is the boost for 7-rays different from 
that of CRs ? The argument is actually very simple. Indeed, for 7-rays, the 
signal is integrated along a line of sight, within a solid angle featured by the angu- 
lar resolution of the telescope. It is therefore very natural to find that the boost 
for 7-rays is a function of the Galactic latitude: at low latitude, one looks toward 
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the GC, where the smooth DM contribution dominates over that of subhalos, and 
the boost is very small; this situation is reversed at high latitude , where the boost 



is much more important, while the absolute flux remains small (Bcrgs trom et al 



1999). For antimatter CRs, one integrates the signal over very different regions 



which are much more concentrated around the Earth, where clumps are not ex- 
pected to dominate. These regions are bounded by the CR propagation scale, 
which depends on both the species and the energy, ft is therefore not surprising 
at all to find results different from 7-rays, even when adopting exactly the same 
subhalo configuration. For comparison, the maximal boost of 20 obtained for CRs 
would give instead a maximum of a few hundreds for 7-rays. Conclusion. We 
have shown that subhalos are not expected to provide an important enhancement 
to the antimatter flux expected from DM annihilation originating from a smooth 
description of the halo. Instead, one of the most important uncertainties affect- 
ing the predictions on the absolute fluxes is the local DM density. We have also 
stressed the origin of the differences between the boosts as computed for CRs and 
for 7-rays: the signals are integrated over quite different regions for those different 
messengers. Finally, we have reminded that the usual description of the Galactic 
DM halo, as very often employed for predictions of indirect astrophysical imprints 
of DM annihilation, is not consistent with the current observational constraints, 
and may lead to over-optimistic expectations. While the effect will be important 
for 7-rays, this has less impact on the predictions for antimatter CRs which are less 
sensitive to the yield from the central regions of the MW. Acknowledgements. 
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Fig. 1. Left: propagation scales for p's and e + 's. The horizontal axis is detected energy 
Ed I injected energy E s for e+'s, while only E / 1 TeV for p's, since the latter do 
not lose energy. Middle: mean boost factors and associated variances obtained for e + 's 
as functions of the energy at the Earth, for different subhalo models. The top curve 
corresponds to the maximal model, which, as expected, is found to increase for decreasing 
propagation scales, but as the variance does. Right: the same for p's, where the energy 
dependence is reversed compared to e + 's. 
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Fig. 2. Left: predictions for the e fluxes for some supersymmetric benchmark models 
as directly calculated from a N-body simulation. Middle: same for p's. Right: baryon- 
subtracted radial velocity curves for the Galaxy, where are also reported the DM contri- 
butions from different density profiles as found in cosmological N-body simulations. 
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